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IDENTIFIABILITY AND GLOBAL STABILITY ANALYSIS ON 
SOME PARTIAL DIFFERENTIAL ALGEBRAIC SYSTEM 


YUSHAN JIANG AND QINGLING ZHANG 


Abstract. We analysis some singular partial differential equations systems(PDAEs) 
with boundary conditions in high dimension bounded domain with sufficiently 
smooth boundary. With the eigenvalue theory of PDE the systems initially 
is formulated as an infinite-dimensional singular systems. The state space 
description of the system is built according to the spectrum structure and 
convergence analysis of the PDAEs. Some global stability results are pro¬ 
vided. The applicability of the proposed approach is evaluated in numerical 
simulations on some wetland conservation system with social behaviour. 


1. Introduction 

Singular systems have abilities in representing dynamical systems in the areas 
of electrical circuits and multibody systems, chemical engineering and economic 
systems, mechanical structure and biological systems. Singular systems are also re¬ 
ferred to as descriptor systems, differential-algebraic systems(DAEs), implicit sys¬ 
tems or generalized state-space systems[?]. A large number of fundamental ideas 
and results based on state-space systems have been successfully extended to sin¬ 
gular systems, such as controllability observability, pole assignment, stability and 
stabilization [?, ?, ?, ?, ?, ?] etc. In recent years, partial differential algebraic 
equations (PDAEs) become an independent field of research which is gaining in 
importance and becoming of increasing interest for both applications and mathe¬ 
matical theory.The PDAEs research areas can be classified into three groups: 

• Index analysis and solvability of PDAEs.(see [?, ?, ?, ?, ?, ?, ?]) 

• Innovative and improved numerical methods to solve PDAEs. (see [?, ?, ?]) 

• Control problems and optimization described by PDAEs. (see[?, ?, ?, ?, ?, 

?]) 

During the last few years there has been a tremendous amount of activity on 
PDAEs. Most practical industrial processes inherently distributed in space and 
time involves the use of PDAEs. Examples include integrated circuits [?], Chemi¬ 
cal Reactor [?] and population systems [?]. Compared with DAEs there exist few 
results on PDAEs reported in literature. 

The structure analysis and solvability on PDAEs can be traced back to [?]. With 
the same method of line (MOL) there existed different spatial indexes between dif¬ 
ferent numerical methods. In order to relate properties of the PDAEs to those 
of the resulting DAEs it is necessary to have a concept of the index of a possibly 
constrained PDAEs. Firstly in [?] the perturbation index was defined on a class of 
liner time invariant PDAEs. On infinite domain the wave solutions of PDAEs was 
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also considered[?]. In[?, ?] there was a systematically analysis about three different 
type of indexes(model, perturbation and algebraic). In [?] a consistent representa¬ 
tion of the solution of an initial boundary value problem for PDAEs was proposed. 
The index involved in the problem is characterized by means of the Fourier and 
Laplace transformations. The index jump was also discussed. The differentiation 
index [?] of some general nonlinear PDAEs on hyper-plane domains is a general¬ 
ization of the differentiation index of DAEs. The differentiation index provide a 
way to determine a Cauchy data on domain surface which must be consistent with 
the PDAEs. For the first time a perturbation index for a singular PDE of mixed 
parabolic-hyperbolic type was computed by [?]. [?] gave some index of PDAEs with 
the sequence of matrices method. It is an effective method for the finite systems. 
But the question is the matrices sequence built in PDAEs are unbounded operators. 
However, it generalizes the Kronecker index in a rather functional analytic manner. 
Other researchers also investigated the index structure of PDAEs systems in the 
area of coupling nonlinear PDAEs systems [?], mixed index [?], index determination 
algorithm [?]. 

Despite the complexities analysis of index on PDAEs, research in the area of 
control problems for PDAEs is relatively scarce. This approach was first applied in 
[?] to design energy based controller of a coupled wave-heat equation systems. In 
[?]stabilizing a coupled PDE-ODE systems with interaction at the interface with 
boundary control was considered. [?] designed a nonlinear saturating control law 
using a Lyapunov function for the averaged model of the switched power converter 
system. Based on the infinite-dimensional state-space representation theory [?] 
addressed the linear quadratic regulator control of the PDAEs. The optimal control 
problem is treated using operator Riccati equation approach. Thought the previous 
methods)?, ?] are derived from DAEs theory. Other optimal control problem)?, ?, ?] 
are considered with parabolic type PDAEs. 

However in these coupled PDE-ODE systems, parabolic distributed PDFs sys¬ 
tems or hyperbolic systems above, the time state variables matrix is reversible and 
the spatial variables performs in one-dimensional interval. Some singular systems 
like parabolic-elliptic partial differential equation can not be direct application of 
the above theoretical study. Our study derived from the article )?] in which the 
search for series solutions of PDAEs with method on DAEs. Motivated by the tech¬ 
nique in )?], we consider some singular partial differential equations systems with 
boundary conditions in high dimension bounded domain. The present work focuses 
on the development of an generalization stability method for a class of PDAEs with 
singular time derivative coefficient matrix. In such a systems the spatial variables 
act on the bounded high dimension domain. 

The organization of the study is as follows: The problem statement for some 
parabolic PDAEs is given in section 2. In section 3 the original PDAEs is described 
as an infinite-dimensional singular systems. Then the state variables expression is 
built with Kronecker-Weierstrass form. And the spectrum analysis is given to show 
the analytical solution of PDAEs is convergence. The dynamical stability properties 
are analysed with Lyapunov method. Finally, in section 4 as an application we build 
some PDAEs model on some wetland conservation system with social behaviours. 
And the stability property of the ecosystem is given by our developed PDAEs 
theory. 
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Notations: K,K" and are the set of real numbers, the n—dimensional 

Euclidean space and the set of all real n x m matrices respectively. || • || denotes 
the Euclidean norm for vectors. For a vector x{t) € K",||a;(t)||oo =: sup ||a:(t)||. 
x{t) G .ifoo if ||a^(i)lloo < 00 . For a symmetric matrix M,M > (<)0 means that it is 
positive (negative) definite. / is the identity matrix.The superscript T is used for 
the transpose. Matrices, if not explicitly stated, are assumed to have compatible 
dimensions. For the convenience, we define the following Hilbert space: 

Jif = {x : ft X [0, +oo) —> K" and ||a :||2 < oo} 
with inner product and o?^-norm respectively defined by 

(a:i,a;2) = {J^xJx2dz},\\x\\2 = . 


2. Description of PDAEs 


We consider the linear partial differential algebra equations system (PDAEs) in 
d—dimensional bounded spatial domain with a state space description of the form 


( 2 . 1 ) 

( 2 . 2 ) 


Ox 

E— = DAx + Ax + Bu, 
ot 


y{t) = [ 

Jn 


Cxdz + V. 


subject to the boundary conditions(BCs) 

dx,{t,z) 


dn 


= 0, z G dfl, L € -dsc 


(2.3) p,x,{t,z) + q,- 
and the initial conditions (ICs) 

(2.4) a;„(0, z) = a:°(z),z G D, It G d/c 


where E,A,D G 


,D > 0, B £ 


X = [xi{t,z) ■ ■ ■ Xn{t,z)\^ G M" is 
the vector of state variables, D C is the spatial domain of the process, z = 
(zi, • • • ,Z(j) G D is the spatial coordinate, t > 0 is the time, A = is the 

Laplace operator, u{t,z) G is the manipulated input vector function, y(t) G 
K"’!' is the measure output, v{t) is the measurement disturbance. Pi,qi, are given 
constants, n is the outward normal vector of x on dfl, dBC C {1, • • • ,n}, die C 
{1, • • • , n}, (z) is the spatial initial state function. 

Remark 2.1. For these systems we show that at least one of the matrices E,D is 
singular. The special case D — 0 leads to DAEs. Therefore in this study we assume 
that D is not a zero matrix. Another special case |if| ^ 0 leads to the PDF-ODE 
systems which some researchers investigated)?]. In contrast to the problems with 
nonsingular matrices E and D, the BCs (2.3) and ICs (2.4) have to fullfill certain 


consistency conditions. Thus, for given PDAEs (2.1) the subsets dsCtdic are de¬ 


termined by consistency conditions (see)?]). In general, the BCs for t ^ dBC and 
the ICs for k ^ djc must be determined with the help of the PDAEs |2.1[ ). With 
different Po,qi values the BCs can be three boundary types: Dirichlet, Neumann 
and Robin boundary conditions. We assume that the BCs are homogeneous for 
simplicity, thought our propose is applicable to the general inhomogeneous bound¬ 
ary conditions. Moreover, we assume that the process in (2.1)-(2.4| evolves on a 
compact set, i.e., {x{t, z),u{t)) £ X x U for all f > 0, where X x U C ffi" x I 
a compact set containing the origin. 


IS 
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3. Decomposition of PDAEs and spectrum analysis 


3.1. Decomposition of PDAEs and infinite singular systems. In this sec¬ 


tion,the PDAEs (2.1)-(2.4) will be rebuilt as a large infinite dimensional systems. 


Consider the linear elliptic operator 

(3.1) Six = DIS.XAx 

in with homogeneous BCs (2.3). We denote a{S) be the spectrum of S that 


is the set of A € C for which (S — XI) is not invertible. More specifically, for the 
operator S, the eigenvalue problem is defined as 

(3.2) 


S(j)j{'Z‘) = Xjfjjji^), (j = 1, • • • ,oo),(/)j(z) € Sf 2 ,n with BCs (2.3). 


denotes the jth eigenvalue and 0j(z) denotes the corresponding 
orthonormal eigenfunction. Obvioursly, the eigenfunctions form an or- 


where Xj S 


thonormal basis for eigenvalue problem (3.2). Furthermore in the next subsection 


we show that a(S) has dscrete spectrum consisting only of real eignevalues with at 
most a finite mumber of positive eigenvalues which is a generalization of [?]. Ap¬ 
plying for PDE theory, we can obtain the following infinite-dimensional differential 
algebric equations systems (IDAEs): 


(3.3) 

(3.4) 

with the initial condition 


EX, = X,X, 

Y, = C,X, + V, 




, 00 


(3.5) 
where 

(3.6) 


Aj(0) = (a;o(z),<;i^-(z)) 


Xj = {x,(l)j),U, = = {v,(l)j{z)),Cj = C [ (l)j{z)dz. 

Jn 


Remark 3.1. One can not guarantee that each system are solvable without consid¬ 
ering the regularity of the matrix pencil set {{sE, Xjl)}^,. Thus in this study we 


assume that for each j the matrix pencil with respect to the system (3.3) and (3.4) 


is regular. Furthermore, all systems have the uniform differential time index)?], 
that is, for all j, the pencil {sE,XjI) is regular and has the Riesz index vj_ j = 
(independent of j). 

3.2. Eigenvalue estimation and property. For the unforced PDAEs 

dx 


(3.7) 


E 


dt 


= DAx -I- Ax. 


with the corresponding IDAEs 

(3.8) BXj = XjXj,j = 1, • • • , oo. 

For every j, firstly let us denote a{E,XjI) be the generalized eigenvalue set of 
systems (3.8) where Xj is defined by (3.2). We also write cr{E) = a{I,E) is the 


usual eigenvalue set. And for every component Xk{t, z) of state vector x, denotes 
the eigenvalue of the Strum-Liouville problem 

(3.9) 


Axk(t,z) + fijXk{t,z) = 0 subject to (2.3). 


The subsequent corollary concerning the eigenvalues estimation of Xj for (3.8) cor¬ 
responding to systems (2.1)-(2.4). Before that, we give two lemmas about the 


spectrum structural property of the scalar elliptic operator. 
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Lemma 3.2 ([?]). With the strongly scalar elliptic operator P (P is strongly elliptic 
in the sense that there exists 7 > 0 such that for all f,, j ^ tICPJ the 

operator defined by = Px + ax = dzj{Pi,j{z)dz.x) + a{z)x in the bounded 

domain G with homogeneous boundary conditions of the form 

dx 

(3.10) — \-P(z)x = 0 on dfl. 

dn 

has discrete spectrum consisting only of real eigenvalues. If a(z) is bounded in It, 
IP can have at most a finite number of positive eigenvalues. 


Lemma 3.3 (The Bauer-Fike Theorem). Assume that A G C”’” is a diagonalizable 
matrix, V G (7"’" is the non-singular eigenvector matrix such that A = VAV~^, 
where A = diag{pi, p, 2 , ■ ■ ■ ,Pn) is a diagonal matrix, Let A be an eigenvalue of 
A-\-T then there exists fi G <j{A) such that |/i — A| < ||k^|| 2 ||T|| 2 ||l^“^|| 2 - 

Corollary 3.4. For the eigenvalue problem the operator P has discrete spec¬ 
trum consisting only of real eigenvalues. And there is a finte number of positive 
eigenvalues,i.e.,if all eigenvalues Xj are ordered that Xj > Aj+i, then there exist a 
finite number g so that Ag+i < 0 and e = < 1 is a small positive number. 


Proof : Considering the Strum-Liouville problems 
(3.11) Axk — p^jXk = 0 with BCs (2.3)(j = 1, • • • , 00 , fc = 1, • • • ,n) 


where denotes the eigenvalue of (3.11). It is true that all pL^ satisfy Lemma 
|3.2[ And from the assumption D > 0 one can get that for every given j if for all 
k = 1, - ■ ■ ,n,p’- < 0 then all eigenvalue of DAx is negative defined. Therefore, it 
is following from Lemma 3.3 that for every j, and Xj satisfy \p^ — Xj \ < c, where 
c is some given constant independed with j. Noticing that -G —oo(j —)• 00 ), the 
conclusion holds. 


Remark 3.5. It should be pointed out that the existing one dimension result [?] 
can not be directly generalized to high dimension case since the mathematical 
complexity property of the Strum-Liouville(L-S) problem on high dimension spatial 
space. Additionally, the zero eigenvaule Xj = 0 yields the trivial constant state 
response x(t,z) = const, which is not interested in our study. 


From Corollarj3.4 there exists some p G such that the operator P with 
respect to (3.2) has discrete spectrum structure 


a{P') — {Ai > X2 > ■ ■ ■ > Xp > Ap_|_i > • • •} 

and Ap+i < 0. Now for the matrix pencil {E, Xjl) with respect to IDAEs ( |3.3[ ) if Sj G 
a(E, Xjl),e G <j{E) then from the definition of eigenvalue the relationship among 
Xj , Sj and e is s^- = ^. Immediately we have the following spectrum structure 
theorem. 


Theorem 3.6. Assume is the spectrum set of the S-L problem (3.2) with 

respect to systems (3.^). e G cr[E) is the eigenvalue of E, then there exist some 
positive number p > 0 such that for every j G ,j > p, the sign of the real part 
of Sj is reversal with e. Consequently, E is a nonnegative matrix is the necessary 
condition for which the matrix pencil family set {{E, Xjl)}^^ is admissible. 
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3.3. State and output response of PDAEs. In what follows we shall assume 


that each pencil {E,XjI) with respect to (3.3) is regular and impulse-free. Intu¬ 


itively, the Kronecker-Weierstrass equivalent form of DAEs[?] can be applied to our 
systems. However, considering the nonnegative diffusion matrix D this equivalent 
transformation leads to an unexpected matrix. And the spectrum analysis theorem 
above is no longer valid. Thus the Jordan transformation is introduced to solve the 
problem. 


For the IDAEs (3.3), there exists non-singular matrix M such that {E, Xjl) can 


be transformed into the following form 


(3.12) 


M-^[E,XjI]M = 


El 

0 


. (A,T) 


where J G is a nilpotent matrix with index no, Ei S R’’^'’ is a nonsingular 

matrix. It follows immediately from the linear singular systems theory)?] that the 
state response of such sytem is 

where Xj^nii{t) are given by 

1 1 

(3.13) Aj,reg(t) = Me^^^i”'*Aj-reg(0)-kM / e^^'®i"'(*"®^Hit7,(s)ds 

Jo 

and 


no-l 


(3.14) 


A,,„ii(t) = -MY, (t). 


i=0 


Hi, i ?2 are the components of M * B corresponding to the partition of the vector 
Xj{t) into the regular and nilpotent parts. Thus the PDAEs (2.1)-(2.4) have the 
state response 


OO 


x{t,z) = } Xj{t)(j)j{z) 


j=i 


and the output response 


y{t) = Y.(c,x,{t) + v,it)) 


j=i 


where Cj,Vj{t), {j G Z+) are defined by (3.6). 


The above state and output reponse discuss provides a generalized systems the¬ 


oretical approach to the PDAEs. From (3.13) we know that if there exists an eigen¬ 


value of El with negative real part then the solution grows exponentially which 
coincides with the ‘explosive solution’ from a mathematical perspective. 

It should be noted that the existing works in[?, ?] rely on the assumptions 
that the derivative matrix E is nonsingular. And for the orthogonal property of 
the eigenfunction set {(()j}^i, some dynamical properties of the PDAEs including 
stability, stabilizability, and detectability can be given through a direct application 
of the LMI technique about the generalized systems theory)?, ?]. The following 
theorem shows the admissible of PDAEs via LMIs. 
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Theorem 3.7 (Admissible Via LMIs). The PDAEs ^3.7\ ) is admissible if and only if 
there exists a matrix P such that the first subsystem of IDAEs (3.S^ is admissible,Le. 

E^P = P^E > 0 , 

Ai(P^ + P) + P'^A + A^P < 0 


where Ai is the maximum eigenvalue of the spatial differential operator (3.1). 


Proof: Noticing that all eigenvalues \j be ordered as \j > A^+i, it can be easily 
deduced from above that 


Xj + P) + P^A + A^P < 0 

holds for every j £ TP'. Then the desired result follows immediately by the admis¬ 
sible theory of the pencil {E,A){see [?] for detail). 


4. Stability analysis 

In view of theorems |3.6| and |3.7| above, the systems stability can be determined 
by spectrum analysis and LMIs. This validity of the above stability analysis relies 
on the convergence properties about the IDAEs. In the following, we propose some 
exponential stability property on the PDAEs by the energy estimation theory. For 
the following considerations, it will be simplest to assume homogeneous Neumann 
boundary conditions. 

Lemma 4.1 ([?]). Let x G Wf (LI),then if p,i is the smallest positive eigenvalue 
of —A on LI (with the appropriate boundary conditions) the following Poincare 
inequalities hold: 

(4.1) ||Vai|p > p,i\\x - x|p, ||Aa:|p > ^i||Va:|p if ^ = Q on 9D; 


(4.2) 

where x = ^ xdz. 


||Vx|p > ^i||x|p if X = 0 on dLl, 


Theorem 4.2. Assume that x(t,z) is a bounded solution of (2.1),{2.2)with ho¬ 
mogeneous Neumann boundary conditions. Assume that E > 0, D > 0, jii is the 
smallest positive eigenvalue of —A on D, di is the smallest positive eigenvalue of 
D and 


g — 2(diAii - mil) ^ Q 


Then 

for a positive constant Ci, and 


\\E\\ 


dx'^ „dx , 

X] -r^E—dz< cie"^* 


i=l 


dz, dzi 


(4.4) [ \\x{t,z) - XM(t)\\dz< C2e 

JQ 

where C 2 is a positive constant, XM{t) = [^ /n z)dz is the spatial average func¬ 
tion, i.e. the state variable vector x(t, z) generated by the PDAEs is exponentially 
stable and asymptotically converge to its spatial average. 
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Proof:By introducing the energy integral(lyapunov function) 

d 

(4.5) ' 

and compute 


/ 5 1 /" ^ dx , 


(4,6) 


= 


dx'^ d'^x 


E 


fdx^ d 

d 

E 


(DAx + Ax) I dz 


2=1 

dx 


/ dx^ ^ dAx dx^ ^ dx 
\ dzi dzi dzi dzi 


dz 


Notice that the Neumann BCs are = 0 on dfl. With the application of divergence 
theorem we have 


(4.7) 


d 

dt 


EL{t) = — [ Ax^DAxdz+ [ 

Jn Jnfri dz, dz. 


The above equality can be estimated as 


(4.8) 


dx^ . dx 


-£;i(<) < -dillAxf + /„ E 7)77^7);:'^^ 


dzi dzi 


where di is the smallest positive eigenvalue of D. According to Lemma 4.1 (4.8) 
implies 


(4.9) 


^i?L(i) < -dimIIV(rf +IIAllIIVxf. 


Noticing that —di/ii + ||A|| < 0, and takeing into account 


dx'^ ^ dx 


(4.10) = 




we obtain that 

(4.11) 

where 


-EUt) < -SEUt) 


^ 2(di/ri - ||A||) ^ 


\\E\\ 


By lemma 4.1 (4.11) implies (4.4). Thus, under the conditions of the theorem, 
spatial oscillations decay exponentially, and the solution asymptotically behaves 
like its spatial average. 
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The sea 


The land 


Wetland area 


Q 


Artificial protection boundary 


5Q 


Figure 1. The idealized spatial domain is a rectangular domain 
with the sea oriented direction zi and coast line direction Z 2 -Tlie 
wetland conservation is closed with no flux boundary conditions 
imposed. 


5. Application to the Coastal Wetland Conservation System with 

SOCIAL BEHAVIOUR 


In this section, we illustrate, through computer simulations, the application of 
the theoretical development given above to some wetland conservation system with 
social behaviour on some plane rectangle domain (see Fig. [^. 

Our model derived from [?] in which a reaction-diffusion system incorporating 
one prey and two competing predator species under homogeneous Neumann bound¬ 
ary conditions was considered. In this model we choose human, birds(predators) 
and their food (prey) as the research objects. The spatiotemporal dynamics between 
predator and their prey with human activity affect in a protected environment can 
be described by the following PDAEs 


(5.1) 


' dxi Xi X2 , . „ 

—- = diAxi + ria;i(l - vw - «! - hiX 3 ),z eil,t> 0, 

ot A'l jN2 

^ = d 2 Ax 2 + r 2 X 2 {-l + ^ 2 ^ - - h 2 X 3 ), z e n, t > 0, 

ot iVi iV 2 

0 = Aa :3 -I- a; 3 , z G fl, 

dxi dx2 dx3 ^ C50 + ^ o 

— = — = — = 0, z e an, t > 0, 

on On On 




where xi{t,z),X 2 {t,z) represent the population of prey (birds) and predator (fish) 
species at time < > 0 and spatial position z G fl respectively; di,d 2 stand for the 
diffusion coefficients of prey and predator species; the prey population follows the 
logistic growth in the absence of predator with the intrinsic growth rate ri and 
the carrying capacity Ni; r 2 is the death rate of the predator with the carrying 
capacity N 2 ]ki, k 2 represent the strength of relative effect of the interaction on the 
two species. 

For a wetland ecosystem, the influence of human can be regard as an invasive 
species and not be affected by other species. Thus the influence of human activi- 
ties(for example,the economic interest) on the two species are represented by hiX 3 
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and h 2 X 3 respectively (in the first two equations of system (5.11) and the human 
population is in a state of free distribution can be described by 


(5.2) z) _ +X 3 (t, z),z e > 0. 

ot 

Since the local human population distribution can reach a time independent dy¬ 
namic balance in a short time, thus the above parabolic equation degenerates to 
the following elliptic equation 


(5.3) 


Axs + X 3 = 0,z G fl. 


Moreover, considering the geographical location affect of sea oriented direction we 
assume that human population X 3 {t, Zi, Z 2 ) = X 3 {t, zi) is independent with the coast 
line direction. And the parameters di, Xi, ki, hi{i = 1,2.) are positive real numbers. 
The systems (5.1) can be rewritten as the following PDAEs 


(5.4) 

where 


Bt 

E— = DAx 
at 


fix) 


E = 



fix) = 


,x = {xi,X 2 jX 3 f' ,D = diag(di,d 2 ,1), 


“ ^ 1 ^ 3 ) 

r2X2i-l + ^ 2 #)- - ^ - h2X3) 
X3 


It is obvious that the positive equilibrium of system (5.1) is the positive solution of 
the following nonlinear equations 


l-|-fc,^-h,.3 = 0,-l. 


^2 ^ - ^ 2 X 3 = 0,X3 = 0. 

1\2 


The positive equilibrium is 


Nijki + 1) A^2(fc2- 1) 
l-|-fcifc2 ’ l-l-fciA:2 


where ^2 > 1 . 

Now we consider the local stable property of the equilibrium. We choose fl := 
{(zi, 22)10 < 2 i < 7 r ,0 < 22 < 1}. To the consistency conditions the initial con¬ 
ditions are taken as ^^(z) = 0 . 3 ,X 2 (z) = 0.3(1 -I- cos( 2 i)). The remaining system 
parameters are chosen as IVi = A ^2 = 1 ti = 2 ,r 2 = 0 . 2 ,/ci = 8,^2 = 18, ,11 = 
diag{2,3,l). By directly computing we have the Jacobi matrix of the linearized 
system at the equilibrium 

/ -18/145 -144/145 -18/ii/145 \ 

Aj = 306/725 -17/725 -17^2/725 . 

Vo 0 1 J 

According to the stability analysis of theorem |4.2| in section 4, we consider the 
local stability property of this PDAEs at the equilibrium in different parameter 
values(see TabQ. Due to the particular choice of the system parameters hi, h 2 the 
state variables xi,a :2 show some different dynamical properties. By increasing the 
value of hi, the value of ndigi — ||A|| can be changed from positive to negative. 
Simulation results for the system in the (t, zi) domain at 22 = 0.5 are depicted 
in Figj^and Fig[^ In Figj^ (a) and (b) the state variables Xi,X 2 spatial-temporal 
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Table 1. Values of the simulation parameters 


Case 

PjII 

difii - ||Aj|| 

II 

II 

O 

1.0069 

1.9931 

hi — 24, /i2 — 0-1 

3.2841 

-0.7159 



Figure 2. Numerical results for the state estimation analysis, (a) 

Time spatial response Xi(t, Zi, 0.5) in the domain; (b) Time 

spatial response X 2 (t, zi,0.5) in the (t,zi) domain; (c)The phase 
trajectory about xi,X 2 {hi = 0.2,/i 2 = 0.1). 

show the corresponding exponentially stable. It shows great difference with the case 
in Fig. (§. Fig.Q shows the different spatial convergence property of the system. 
Ecologically the stable prey-predator relationship stands when the diffusion rates of 
prey and predator diag{di,d 2 ,1 ) are very high with low human affect rates hi,h 2 
which correspond with the experience. 


Conclusion 

In this study we have studied the problem of some paraboli-elliptic type PDAEs 
in high dimensional domain. With the decomposition ideas derived from PDE 
theory we built the IDAEs to reconstruction of PDAEs. Some spectrum theory 
result of the IDAEs corresponding the PDAEs are proposed. A exponential stable 
result on the PDAEs is presented through the energy estimation about the state 
variables under the homogenous Neumann boundary conditions for the positive 
diffusion matrix. Finally as an application we built some wetland conservation 
model with social behaviour. The numerical results show the effectiveness of this 
development. 
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Figure 3. Numerical results for the state estimation analysis.(a) 
Time spatial response Zi, 0.5) in the it,Zi) domain; (b) Time 
spatial response X 2 (t, Zi,0.5) in the {t,zi) domain; (c)The phase 
trajectory about xi,X 2 {hi = 24 ,/12 = 0.1). 



Figure 4. Numerical results on spatial convergence of xi,a;2. 
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